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Abstract. Generation of the large-scale coherent vortical structure in homogeneous shear 
flow couples dynamical processes of energy and enstrophy production. In the large rate of 
strain limit, the simple estimates of the contributions to the energy and enstrophy equa- 
tions result in a dynamical system, describing experimentally and numerically observed self- 
sustained non-linear oscillations of energy and enstrophy. It is shown that the period of these 
oscilaltions is independent upon the box size and the energy and enstrophy fluctuations are 
strongly correlated. 
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Due to its seeming simplicity, the problem of homogeneous shear flow has widely been 
used as a benchmark for nu merical and experimental tests of various closures for turbulence 
modelling. All early closures were based on the Kolmogorov ideas developed for statistically 
steady isotropic and homogeneous small-scale turbulence interacting with the non-universal 
large-scale flow-field. It became clear that to validate this physically appealing concept, one 
had to verify and understand the symmetries and other statistical properties of the small- 
scale velocity fluctuations in the real-life flows. This was the main focus of the experimental 
studies of homogenious shear flow [l]-[6]. The interest in this model flow is also related 
to the recent numerical investigations which revealed coherent structures resembling those, 
responsible for turbulence production, in the wall- sheared flows [7]- [13]. This system has 
also often been used for calibration of various constants in semi-empirical turbulence models 



The problem is formulated as follows: consider a flow in a cube of a side a, so that 
—a<Xi< a. The velocity field 



with the imposed mean velocity < v >= U(y)ei = Syei (the definition of the averaging 
operation will be introduced below). The vorticity is defined then: 



The equations of motion for the fluctuating components of velocity and vorticity are 
(density p = 1): 



[15],[16]. 




(1) 



f2 = —Se 3 + uj 



(2) 



d t u + u ■ Vu = -Sve 1 - Vp - U(y)d x u + z/V 2 u, 



(3) 



d t u) + u • Voj = uj • Vu — S(—d z we 3 + d x we\ + d z v e 2 ) — Syd x u 



(4) 



and 



V-u = 



(5) 
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The x, y and ^-components of velocity field are denoted hereafter as u, v and w, respectively. 
Let us define the averaging operations: 



F =< F(x, t) >= -- / / rfxrftF(x, t) 
1 V Jo Jv 

in the limit V = a 3 — > oo; T — > oo. The statistically steady state is assumed here. The 
spacial averaging is defined as 

F(t) = Fj^t} = ^ |dxF(x,f) 
The kinetic energy equation is {y — > 0): 

d t /C + \uiViU • = -r^S" - VjpMl - £^ (6) 
I 

with r ra = ml. The contribution U (y) = due to the symmetry of the problem. 

Since in a homogeneous flow all spacial derivatives of the mean properties are equal to 
zero, the modelling is reduced to investigation of the time evolution of turbulent kinetic 
energy K, — uf/2 and dissipation rate S = v(diUj) 2 . 

For small perturbations from isotropic and homogeneous state (S — > and pu = pv = 0), 
the typical turbulence models, based on the equilibrium ideas are [15], [16]: 

dtK = -TijSij - S (7) 

and 

£ £ 2 

d t £ = -C el TijSij— - C e2 -^r (8) 

where the Reynolds stress = ujv]. The coefficients C ei = 0(1). A simple expression, valid 
at the long times t > K/£, 

T i:j ~ -v T Sij (9) 

with turbulent viscosity ut oc K 2 /£, closes the set of equations (7)-(9) and defines the so 
called K — £ model, widely used in engineering for modelling the not-too strongly sheared 
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flows. The unknown magnitudes of the proportionality coefficients are typically determined 
in a following way. Consider a flow with 5 = 0. The unknown coefficient C e2 can be found 
from comparing the analytic solution of the simple equations (7)-(9) with experimental and 
numerical data. The same flow can also be used to test the results of analytic theories 
[15]-[16]. If S 7^ 0, solution of (7)-(9) is not easy and the coefficient C el can be found from 
comparison with the data. The solution of equation (7)- (9) with the fixed values of the 
coefficients showed a close to exponential long-time growth of turbulent kinetic energy in a 
good agreement with the outcome of direct numerical simulations [15]- [16]. If the shear is 
imposed on a decaying isotropic turbulence at t — 0, the observed [16] initial, short- time 
decay of kinetic energy is readily explained by the fact that turbulent viscosity 



which is small at short times. 

Recent numerical experiments revealed a much more complicated picture. Driven by a 
very strong shear (the criterion is derived below), in the long-time limit, the system developed 
a limit cycle -like strong fluctuations of the total kinetic energy about the mean value < K > 
[7]- [13]. The amplitude of these fluctuations was up to two-three times that of < fC >. Similar 
effect was observed by Borue et. al. [14] in a three-dimensional Kolmogorov flow driven by 
a steady forcing f = (0, 0, cos(x)). Elucidation of the origin of these oscillations is the goal 
of this paper. 

The physical process observed in both homogeneous shear and Kolmogorov flows can be 
described in two steps: first, the shear generates both kinetic energy and vortical structures 
leading to the access of the energy production. Then, the structures become unstable and 
rapidly disappear with the energy dissipation taking over. The process repeats itself. The 
evolution of kinetic energy and enstrophy fluctuations in 3D Kolmogorov flow, conducted 
by Borue et al [14], revealed extremely strong correlation: the sharp spikes in the enstrophy 
and energy time-signals were almost simultanious with a slight time-lag, thus suporting the 
importance of coherent vortical structures in the process. 

At the long times the numerical homogeneous shear flow problem (l)-(3) has two very 
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important features. We can see from the equation of motion that the flow, defined on a cube, 
cannot be periodic in space. Second, the integral scale C in this situation is not a dynamic 
variable which is a function of /C and £, but prescribed by the box size, so that C ~ a. This 
puts strong constrains on the modelling of various contributions to the equations (l)-(3). 

Now, we would like to establish the main characteristic length-scales. The non-universal 
velocity fluctuations belong to the range of scales a ps C < r < r c with the cross-over scale 
r c ks \[EjS^ are dominated by powerful anisortopic coherent structures (vortices). The 
univesal range, populated by the more or less isotropic excitations, spreads over the interval 
r c < r < rd ~ (v 3 /£)^. The Kolmogorov spectrum can be expected in the range with the 
total energy of quasy-isotropic fluctuations 



J r d 

We can see that the inertial range shrinks to zero when the strain rate becomes large. 
This fact, noticed in Ref. [13], defines the strong shear regime. In strongly anisortopic flow, 
the simple expression (9), is invalid. 

First , let us consider the equation for ui 2 : 



Due to powerful, shear-generated coherent vortical structures, oj y uj x = 0{uj 2 ) and the con- 
tribution involving longitudinal derivative d z w can be neglected. The simple dimensional 
considerations lead to: 




(10) 





u> ■ ui ■ Vu = 0{u rms u) 2 1 £); S{uj x d x w + uj y d z v) = 0(Su 2 ) 



and denoting A(t) = u 2 / S 2 we have: 




(12) 



with all coefficients 7, a and a = 0(1). 
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To estimate r uv , we observe that in the limit of interest (see below) the only relevant 
time scale is r c ~ 1/S. From the equation (3) we have an estimate for the stress: 



T uv = u(t)v(t) « u(t) f d\(-u(A) • Vu(A) - V„p(A)) (13) 

where the "initial condition r u „(i — r c ) was neglected for simplicity (see below). The dimen- 
sional estimate gives: 

Tuv = m sc (14) 

where B(t) is an "anisotropy factor" or "order parameter" characterizing the varying in time 
strength of the coherent vortical structures. The appearence of this factor is natural (see 
below) since in an isotropic, non-sheared flow lacking coherent vortices B(t) = const = 0, 
while in the strongly sheared flow these structures contribute to the energy production. The 
dissipation rate in shear flows is estimated as £ = blC? / C with the coefficient b ~ 1 leading 
to a model equation: 

K t (t) * (B(t) - b) 3 ^- (15) 

or, introducing y = \/lC: 

2d t y(t) = ^^y 2 (16) 

By the virtue of (2), the mean value of vorticity in homogeneous shear flow is < Q >= S. 
Thus, the natural measure of the strength of the strongly anisotropic fluctuating coherent 
vortical structures is the ratio uJ 2 / 'S 2 . Based on these considerations , we set A(t) = B{t). 
This result can be derived in the third-order of the iteration procedure of the expression (13). 
Indeed, inserting an unknown initial condition into (13), one can use r uv (t — r c ) = r as a zero 
order solution. Then, after simple resummation (neglecting the first -order contributions) 
we obtain: 



T U v{t) = r + f u(t)u(\) ■ Vv(\)+ IN u(t)u(\) ■ Vu(A') • Vu(A") • Vv(\")d\d\'d\ 

Jt — T c J J J 



This expression immediately gives (14) with B(oj) « provided the derivatives diUj ~ 
uj. This approximate derivation is given here to demonstrate the mechanism of vorticity 
appearence in the expression for the reynolds stress r uv . It will become clear below that the 
power of vorticity in the expression B{oj) oc uj n is unimportant. 

The fact that the anisotropic ordered structures can influence the magnitude and even 
sign of the energy fluxes is known for a long time. The rigorous linear stability analysis 
developed in [17]- [18], showed that even in three-dimensional flows the strongly anisotropic 
structures (basic flows), are capable of reversing the sign of the energy flux, due to the 
"negative viscosity" effects and lead to substantial growth of a small large-scale perturbation. 
In the opposite limit of the isotropic basic flows, the theory showed generation of positive 
effective viscosity and acceleration of the energy dissipation. These features are incorporated 
in the model (12), (15), (16): indeed, we see that if A(t) > b, the energy grows, while, when 
A < b, it decays. The model equation (12) includes the well-known process of the vortex 
break -down: when v 2 » u 2 £ 2 , the instability leads to the vortex disappearence. 
Defining dimensionless variables Z = jfe, Z = ^ and T = St gives: 

2d T Z{T) = a(A(T) - b)Z 2 {T) (17) 

and 

d T A{T) = -i(Z(T) - Z )A(T) (18) 

where Z > = const, related to the mean amplitude of Z(t). 

The equations (17), (18) have a steady-state solution Z(T) = Zq and A(T) = b. A simple 
linear stability analysis shows periodic solution when the amplitude of the perturbation is 
very small. The numerical solutions of quations (17), (18), presented on Figs. 1-6, revealed 
strong non-linear oscillations. All calculations were perforemed with Mathematical™ . In 
a wide range of parameter variation, the system generates non-linear oscillations with the 
shape depending upon initial values Z(0) and A(0). For a given set of parameters the 
frequency of oscillations is proportional to the strain rate S. 
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For the initial values of Z(0) and A(0) ~ 1, the solution shows reasonably smooth 
oscillations with Z and A being somewhat out of phase (see Figs. 1,2). The result supports 
a general physical picture of the anisotropy A(t) (order parameter) and energy growing 
(decaying) together with some time-lag. The energy fluctuations are by a factor 2-3 larger 
than S 2 C 2 . When the initial energy was doubled to Z(0) = 3., the oscilaltions became much 
less symmetric with the steeper energy grows (Figs. 3,4). The crucial role of the "order 
parameter" A(0) is demonstrated on Figs. 5,6 corresponding to Z(0) = 2 and A(0) = 0.1. 
We can see the the formation steep shock-like structures, somewhat resembling turbulence- 
production bursts. 

In the range of large Z(0) and very small A(0) << 1, the solution blows up, indicating 
the unphysicality of these initial conditions corresponding to the large energy fluctuations 
(u rms (0) » SC) and small anisortopy (order) parameter {u 2 « S 2 ). 

To discuss the above results, let us look at this work from a somewhat different angle. The 
Kolmogorov relation S 3j0 (r) = (u(x + r) — u(x)) 3 oc r, is a statement about constantcy of the 
energy flux for inertial range wave numbers fc~l/r>>l/£of isotropic and homogeneous 
turbulence. As r — > C, the structure function S^o(r) — > 0. In strongly anisotropic flows 
with the integral scalle C ~ a, this is not so: depending on the spacial distribution of 
velocity (vorticity), the moment S%{£) ^ 0. If vorticity (enstrophy) is an "order" parameter, 
characterizing deviations from isotropy, then S${C) ~ B(lu)K,z / C where B(u) — > when the 
strength of the structures diminishes. This qualitative statement is supported by the well- 
known fact that the velocity field, generated by the vortex v(r) oc F<f>(r), where the circulation 
T = OijjjC 2 ). Combined with the equation for the enstrophy, the two relations (17), (18) form 
a dynamical system leading to strong fluctuations of both energy and enstrophy. The shape 
of the function B(uj) does not seem to influence the qualitative aspects of the process: the 
model (17), (18) is invariant under transformation A(t) — > A n with a simple rescaling of time. 

All this is valid when C ~ a. If this is not so, the magnitude of the fluctuations must 
substantially decrease. Indeed, if a >> C, then we are dealing with N = (^) d independent 
systems. Here d is the force dimensionalty. Since the phases are crucially important , we 
expect the amplitude of the fluctuations to decrease as 1/a//V. This can easily be tested on an 
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example of 3D Kolmogorov flow in a box with the side a driven by the force f = (0, 0, cos(f )) 
by varying the forcing scale. 

To conclude: based on the equations of motion and some physical considerations, we 
propose a dynamic model, coupling vorticity (enstrophy) and energy fluctuations in a ho- 
mogeneous shear flow. This model generates strongly correlated self-sustained oscillations 
of both enstrophy and energy similar to those observed in eperiments and direct numerical 
simulations. The calculated time -lag is similar to that observed in a numerical study of 3D 
Kolomogorov flow by Borue et al [14]. 

It is not yet clear if, properly parametrized, this simple model can mimic turbulent bursts 
which are at the core of the energy production in turbulent wall flows. In case of a positive 
answer, the model of this kind can serve as a boundary condition ( "wall function" ) for turbu- 
lence simulations, neglecting the detailed consideration of dynamics of the viscous sublayer. 
The achieved computational economy makes this aspect of the work worth persuing, 
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Figure 4: Parametric plot Z 2 (T) (horizontal) vs A(T) 
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Figure 5: Time -evolution of Z 2 (T) oc K,(T) (higher amplitude curve) and A(T) vs T. Z Q 
A(0) = 0.1 
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Figure 6: Parametric plot Z 2 (T) (horizontal) vs A(T) 
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